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Abstract 

We propose a new method to consider D-brane probes in thermal backgrounds. The method 
builds on the recently developed blackfold approach to higher-dimensional black holes. While 
D-brane probes in zero-temperature backgrounds are well-described by the Dirac-Born-Infeld 
(DBI) action, this method addresses how to probe thermal backgrounds. A particularly im- 
portant feature is that the probe is in thermal equilibrium with the background. We apply 
our new method to study the thermal generalization of the Blon solution of the DBI action. 
The Blon solution is a configuration in flat space of a D-brane and a parallel anti-D-brane 
connected by a wormhole with F-string charge. In our thermal generalization, we put this 
configuration in hot flat space. We find that the finite temperature system behaves qualita- 
tively different than its zero-temperature counterpart. In particular, for a given separation 
between the D-brane and anti-D-brane, while at zero temperature there are two phases, at 
finite temperature there are either one or three phases available. 
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1 Introduction 

The Dirac-Born-Infeld (DBI) action is a low energy effective action for D-brane dynamics 
obtained by integrating out the massive open string degrees of freedom pQ . The first example 
in string theory where the full non-linear dynamics of the DBI action was exploited is that of 
the Blon solution [21 [3] . In the Blon solution the configuration in which an F-string ends on 
a point of a D-brane is resolved into a smooth solution of the DBI action where the F-string 
dissolves into the D-brane. The D-brane carries a world-volume electric flux which is the 
F-string charge. The new phenomena at play are that one can describe multiple coincident 
F-strings in terms of D-branes and furthermore that the F-strings go from being a one- 
dimensional object of zero thickness to be "blown up" to a higher-dimensional brane wrapped 
on a sphere. These phenomena are captured thanks to the non-linear nature of the DBI action. 
Based on these phenomena, many important applications of the DBI action were found in the 
context of the AdS/CFT correspondence For gravitons satisfying a BPS bound by moving 
on the equator of the of the AdSs x background it was found that they blow up to become 
Giant Gravitons, D3-branes wrapped on three-spheres, for sufficiently large energies [HEIIT]. 
Another interesting application is the Wilson loop, originally considered in the AdS/CFT 
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correspondence using the Nambu-Goto F-string action |^i9j. Here, the "blown up" version 
for a Wilson loop in a high-dimensional representation has been considered using the DBI 
action, either for the symmetric representation using a D3-brane [10] or the antisymmetric 
representation using a D5-brane [11]. 

The success of using the DBI action to describe D-branes probing zero-temperature back- 
grounds of string theory motivated the application of the DBI action as a probe of thermal 
backgrounds, particularly in the context of the AdS/CFT correspondence with either ther- 
mal AdS space or a black hole in AdS as the background Ll2j. Applications include meson 
spectroscopy at finite temperature, the melting phase transition of mesons and other types of 
phase transitions in gauge theories with fundamental matter (see |13j for early works on this) . 
Furthermore, the thermal generalizations of the Wilson loop, the Wilson-Polyakov loop, in 
high-dimensional representations were considered [14^ I15j 

However, the validity of using the DBI action as a probe of thermal backgrounds is not 
clear. In general the equations of motion (EOMs) for any probe brane can be written as 
[Ml [IT] 



K^.PT"^ = J . pP (1.1) 



where Tab is the world-volume energy-momentum (EM) tensor for the brane, Kab^ is the 
extrinsic curvature given by the embedding geometry of the brane and the right hand side, 
J ■ FP, represents possible external forces arising from having a charged brane that couples 
to an external field. In the applications of the DBI action as a probe of thermal backgrounds 
the D-brane is treated as if the temperature of the background does not affect the physics 
on the brane. Therefore, the EM tensor that enters in the EOMs (1.1) is the same as in the 
zero-temperature case. However, there are degrees of freedom (DOFs) living on the brane 
that are "warmed up" by the temperature of the thermal background, just like if one puts a 
cold finger in a big bathtub with warm water. The thermal background should thus act as a 
heat bath for the D-brane probe and the system should attain thermal equilibrium with the 
D-brane probe gaining the same temperature as the background. Because of the DOFs living 
on the brane this will change the EM tensor of the brane and thus in turn change the EOMs 
(1.1) that one should solve for the probe brane. 

In this paper we study the thermal generalization of the Blon solution. This serves as 
a test case to study D-branes as probes of thermal backgrounds. The Blon solution is a 
solution of the DBI action for a D-brane probing ten-dimensional flat space-time, the D-brane 
world-volume having an electric flux interpreted in the bulk as an F-string. We shall instead 
consider ten-dimensional hot flat space as our background. The challenge is that one does 
not know what replaces the DBI action, which is a low energy effective action for a single 
extremal D-brane at weak string coupling, when turning on the temperature. However, in 
the regime of a large number of coinciding D-branes we have an effective description of 
the D-branes in terms of a supergravity solution in the bulk when QgN ^ 1. Using this 
supergravity description one can determine the EM tensor for the D-brane in the regime of 
large A'^. This EM tensor will then enable one to write down the EOMs (1.1) for a non- 
extremal D-brane probe in the regime of large A^. This results in a new method in which one 



^Unlike in the D5-brane case, it seems that for the D3-brane case, which corresponds to a totally symmetric 
representation of the Polyakov loop, there is no solution [IHIT5] . 
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can replace the DBI action, which is a good description of a single D-brane probing a zero- 
temperature background, by another description that can describe coincident non-extremal 
D-branes probing a thermal background such that the probe is in thermal equilibrium with 
the background. 

The idea behind the new method for non-extremal D-branes probing thermal backgrounds 
comes from the so-called blackfold approach [TTl [181 HH] which recently has been developed 
in the study of higher-dimensional black holes. Indeed, the conceptual and technical devel- 
opments of the blackfold approach to black holes provide the building blocks that we need in 
this paper. The blackfold approach provides a general description of black holes in a regime 
in which the black hole approximately is like a black brane curved along a submanifoJd of a 
background space-time (hence the name "blackfold"). This regime entails in particular that 
the thickness of the black brane is much smaller than the length scale of the embedding ge- 
ometry]^ The general idea of the blackfold approach is thus that one integrates out the brane 
thickness scale and obtains an effective description in terms of the large scales over which the 
brane varies, such as the scale of the embedding geometry. Doing this one obtains an effective 
description with the world- volume EM tensor of the brane being that of a fluid and with the 
dynamical principle being the conservation of the EM tensor. The EOMs resulting from this 
consist of the conservation of the EM tensor on the world-volume along with the extrinsic 
EOMs for the blackfold of the form (1.1). To leading order the brane can be regarded as a 
probe brane that does not backreact on the background geometry]^ Thus, as we shall see, one 
can parallel the probe approximation in the blackfold approach with the probe approximation 
that the DBI action assumes and the only difference in the extrinsic EOMs (1.1) is that one 
should replace the DBI EM tensor with that of the fluid EM tensor for the black brane. 

Before turning to finding the thermal generalization of the Blon solution, we begin this 
paper in Section [2] by reviewing the original Blon solution. As part of this we introduce the 
physical set up of the Blon such as how we embed the Blon solution in flat ten-dimensional 
space-time. We restrict ourselves in this paper to consider D3-branes. We review that the 
Blon solution either takes the form of an infinite spike with an F-string charge coming out 
of the D3-brane, describing a number of coincident F-strings ending on the D3-brane, or the 
form of a D3-brane and a parallel anti-D3-brane connected by a "wormhole" with F-string 
charge. The latter solution would correspond in the linearized regime to a string stretched 
between a D3-brane and an anti-D3-brane. 

In Section[3]we consider the question of D-brane probes in thermal backgrounds in general. 
We approach this in Section |3.1| by first reconsidering the general form of the EOMs for the 



DBI action and show that they can be cast in the form (1.1) with the world- volume EM 
tensor being given from the DBI action. Since (1.1) are the general EOMs for a probe brane 
all we need to do in order to obtain the EOMs for a thermal D-brane probe is to find the 



^The brane thickness scale is given by the length scale over which the brane backreacts on the surrounding 
space-time. 

''obviously, if the distance to the brane is of order the brane thickness scale one would observe a backreaction 
to the surrounding geometry from the brane. However, we are precisely integrating out the brane thickness 
scale, thus the probe approximation assumes in this case that the backreaction from the brane is sufficiently 
small to be neglible for the large scale physics described by the extrinsic EOMs I.e. only when considering 

corrections to the probe approximation one should start taking into account the backreaction of the brane on 
the background geometry when computing the extrinsic curvature tensor Kab'' ■ 
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corresponding EM tensor. This we attack in Section 3.2 Since the Blon solution includes an 
F-string charge on the brane we find the EM tensor for the non-extremal D3-F1 brane bound 
state. We subsequently explain how this can be used to write down the EOMs for a thermal 
D-brane probe by employing the blackfold approach. We find in particular that the EOMs for 
the thermal D-brane probe in the extremal limit are the same as for the DBI action probe, 
if one sets the number of D3-branes to 1, while keeping the ratio of the F-string charge 
relative to the D3-brane charge fixed. Finally, in Section 3.3 we give a detailed argument why 
a D-brane probe in a thermal background is not accurately described by the DBI action. 

In Section [4] we consider the Blon solution in the background of hot flat space. We 
first write down the EOMs (1.1) explicitly for the setup for the Blon solution introduced 
in Section [2] This results in a single equation for the profile of the solution. Based on the 
blackfold approach, we propose an action from which the EOMs can be derived. For our setup 
this gives the same equation as before, but written in a form that makes it readily solvable. 
Using this we find the general solution for the D3-F1 blackfold configuration that generalizes 
the Blon solution of the DBI action to non-zero temperature. 

In Section [5] we examine more closely the brane-antibrane wormhole solution at finite 
temperature. This solution consists of N D3-branes and N anti-D3-branes, separated by a 
distance A, connected by an F-string charged wormhole with minimal sphere radius aQ. We 
consider A as a function of fio when ao is varied for a given temperature and compare in this 
way our configuration at non-zero temperature with the corresponding extremal configuration. 
At zero temperature, the wormhole solution is characterized by a thin branch with small ao 
and a thick branch with large ctq for fixed A. We discover that, when the temperature is 
turned on, the separation distance A between the brane-antibrane system develops a local 
maximum in the region corresponding to the zero temperature thin branch. This is a new 
feature compared to the zero temperature case. The existence of this maximum gives rise to 
three possible phases with different o"o for a given A. For small temperatures and/or large ao 
the A as a function of fio resembles increasingly closely the zero temperature counterpart. 

We end in Section [6] with the conclusions and discussions of further directions. 

In a forthcoming paper |2Uj we investigate the thermodynamics of our newly found solution 
by analyzing the physics of the three different phases found in Section[5]in the brane-antibrane 
wormhole configuration and by examining whether it is possible to construct a thermal gen- 
eralization of the infinite spike solution. 



2 Blon solution of DBI revisited 

At low energies, the effective theory describing D-brane dynamics is given by the DBI action. 
The DBI theory is a non-linear theory of a U{1) gauge field living on the world- volume and a 
set of neutral scalars describing transverse fiuctuations. In the linearized regime, this reduces 
in particular to Maxwell electrodynamics for the gauge field while the scalars are free and 
massless. In this linear approximation there exist world- volume solutions describing Maxwell 
point charges with delta function sources. In the bulk theory these are interpreted as F-strings 
ending on a point charge on the D-brane. However, this picture is changed when one takes into 
account the non-linearities of the DBI action. It has been found that the brane curves before 
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one reaches the point charge. The solution of DBI theory describing this is called a Blon, and 
its bulk interpretation is that of F-strings dissolving into the D-brane j2l [3] . Furthermore, 
these Blons solutions can describe a configuration of two parallel D-branes (with one being 
an anti-brane) connected by F-strings. 

In this section we review the Blon solution since our aim in this paper is to describe 
its thermal generalization. This will also enable us to introduce the setup and provide the 
results to which we can compare the thermal solution in the zero temperature limit. We will 
first briefly review how the solution is obtained from the DBI action using the Hamiltonian 
approach and subsequently describe the corresponding spike and wormhole solutions. 

We consider in this paper only the D3-brane Blon solution, and its thermal generalization. 
However, all our results and considerations can readily be extended to general Dp-branes. 

2.1 DBI action and setup 

Consider a D3-brane embedded in a 10- dimensional space-time with the only background flux 
turned on being the Ramond-Ramond five-form flux. We furthermore assume the background 
to have a constant dilaton. The DBI action for the D3-brane then takes the form 

/dbi = -Td3/ dVV- det(7,b + 2TTp,Fab) + Tds / P[C(4)] (2.1) 

J w.v. J w.v. 

where the integrals are performed over the four-dimensional world- volume. Here we have 
defined the induced world- volume metric 

lab = g^udaX^'da" (2.2) 

where g^u is the background metric, X^(cj") is the embedding of the brane in the background 
with cr" being the world- volume coordinates, a, 6 = 0,1,2,3 are world- volume indices and 
/u,z/ = 0,1,..., 9 are target space indices. Furthermore, Fab is the two-form field strength 
living on the D-brane, C(4) is the RR-four form gauge field of the background and P[C(4)] is 
its pull-back to the world-volume. Finally, the D3-brane tension is Tds = [{2iTY'gslf[~^ where 
gs is the string coupling and Ig is the string length. 

Embedding 

To describe the Blon we specialize to an embedding of the D3-brane world volume in lOD 
Minkowski space-time with metric 

6 

ds^ = -dt^ + dr^ + {de"^ + sin^ Odcj)'^) + ^ dx^ (2.3) 

i=l 

without background fiuxes. Choosing the world volume coordinates of the D3-brane as 
{cT°, a = ... 3} and defining r = cr'^, o" = cr^, the embedding of the three-brane is given 
by 

t(a'^)=r, r(cT") = f7, xiia") = z{a) , eia"") = a\ (/){a'') = a\ (2.4) 

and the remaining coordinates Xj=2..6 are constant. There is thus one non-trivial embedding 
function z{a) that describes the bending of the brane. The induced metric on the brane is 
then 

-iabda^da^ = -dr^ + {} + ^'i'^f) ^^'^ + (^^^ + sin^ ^dcj)'^) ■ (2-5) 
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so that the spatial volume element is (iV(3) = y^l + z' {a)'^a^ (1^1,(^2) ■ 

To get the appropriate F-string flux on the brane we turn on the world- volume gauge field 



strength component Fqi. With this, the DBI action (2.1) gives the following Lagrangian 
L 



(27r/2Foi) 



2.6) 



Note that we assumed Fqi to depend only on a since this is required for spherical symmetry. 
Boundary conditions 

We have two boundary conditions on the Blon solution. The first one is 

z{a) — )• for a 



oo . 



(2.7) 



This condition ensures that far away from the center at r = the D3-brane is flat and infinitely 
extended with xi = 0. Decreasing a from oo the brane has a non-trivial profile xi = z{a). In 
general we have a minimal sphere with radius ctq in the configuration. For a Blon geometry 
z{a) is naturally a decreasing function of a where at a = ao the function z(a) reaches its 
maximum. Thus, a takes values in the range from ctq to oo. At ao we impose a Neumann 
boundary condition 



z'(a) —7- — oo for a 



The rationale of this condition is that if z{aQ) < oo the brane system cannot end at (r = 
(To,xi = z{ao)) because of charge conservation and this boundary condition, as we describe 
below, enables us to attach a mirror of the solution, reflected in the hyperplane xi = z{ao). 
In line with this, we define 

A = 2z(cJo) (2.9) 

as the separation distance between the brane and its mirror. Fig. [T] illustrates our setup and 
the definitions of o"o and A. 

20-0 




Figure 1: Illustration of the setup, showing the embedding function z{a) and the definition of the 
parameters ctq and A. 



2.2 Blon solution 



We now consider the Hamiltonian corresponding to the Lagrangian (2.6). To derive this we 
need the canonical momentum density 4:'Ka'^JI{a) = 6L/5{drAi) associated with the world- 
volume gauge field component Ai. Using Fqi = drAi — O^Aq this gives 

Uja) = _Ji <^)'^°^ (2.10) 
l + ^''-(27rZ2Foi)2 
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so that the Hamiltonian can be easily constructed as 

Hbbi = daa^n{a)drAi{a,T) - L = An J da [a'^U{a)Foi - d„{a^U{a))Ao{a)] - L 

(2.11) 

Here we have in the second step integrated by parts the term proportional to daAo, showing 
that Aq can be considered as a Lagrange multiplier imposing the constraint So- ((7^11(0")) = 
on the canonical momentum. Solving this constraint gives 

k Tdsk 



(2.12) 



where k is an integer, Tpi = l/27rZ^ is the tension of a fundamental string and we have defined 
K = = kTrgJg. Using ( 2.12| ) in (2.11) the Hamiltonian becomes 



Fdbi = 47rrD3 J da^/l + z'{ayFBBi{(r) , Fdbi(ct) = 1 + ^ 



(2.13) 



The resulting EOM for z{(7), obtained by varying (2.13), is 

[,/1 + z'iayJ 

Solving for z{a) subject to the boundary conditions stated above, we obtain 



(2.14) 



. f i^DBl(^)' ^ 



(2.15) 



where we recall that do is the minimum value of the two-sphere radius a. The explicit solution 



for z can be obtained by integrating the expression in (2.15) 



z a 



- a, 



(2.16) 







In particular for non-zero ao this represents a solution with a finite size throat, as illustrated 
in Fig. [2} 



Figure 2: Sketch of solution with finite size throat. 



From the expression in (2.15) we can compute the energy density by evaluating the inte- 



grand of (2.13), yielding 

dH _ 
da 



47r Td3 ^(1 + ^'(ct)2) (ct4 + k2) = 47r Tds 



a^ + K^ 



cr-*-a?, 



(2.17) 



From this, dividing by the derivative of the solution z'{a), we compute the energy density 
along the brane 



dH 



1 dH 



dz z'{a) da 
which is finite for a in the range [do, oo). 

Spike solution 



47rTD3 



(2.18) 



For do = the integral in Eq. (2.16) gives the Coulomb-charge type of solution 

z{a) 



a 



(2.19) 



i.e. the spike solution (see Fig. [s]). In Ref. [2] it was shown that the energy corresponding 
to this solution is the energy of a fundamental string of a given length. This was done by 



comparing the integral appearing in (2.13) with the explicit form of the solution z{a) at a 



point near the end of the spike, in the linear approximation and with a suitable regularization 
of the integral providing the energy. In the non-linear case, and also to avoid divergences, it 



is more convenient to compute the energy density along the brane as in (2.18). In particular 
setting cj = (To = in (2.18) we find that the energy density at the tip of the spike is given by 



dH 

dz 



(2.20) 



=(T0=0 



where we used k defined below (2.12). We thus find that this is the tension of the fundamental 



string times the number of strings k, as expected. 



J 



Figure 3: Sketch of the spike configuration. 



Wormhole solution 



For (To = we showed above that the solution (2.16) corresponds to a spike. However, as 



explained in Sec. 2.1 for more general values of ao one can use the solution to construct a 
configuration representing strings going between branes and anti-branes [2|, to which we refer 
as the wormhole configuration (see Fig. |4|. 

Obviously, a system of a D-brane separated from an anti-D-brane is unstable since the 
branes attract each other both gravitationally and electrically. However, the time scale of 
this is very large for small string coupling gg <^ 1 since the tension of a D-brane goes like 
l/ffs while the gravitational coupling goes like gg. Indeed, for a Dp-brane and anti-Dp-brane 
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Figure 4: Attaching a mirror solution to construct a wormhole configuration. 



system the time scale t is of order 



(2.21) 



where A is the distance between the brane and anti-brane. 



Turning back to the wormhole configuration, the separation (2.9) between the D3-brane 



and anti-D3-brane can be computed from (2.16) and is given by 

2Vvrr(|)74T^ 



2z(ao) 



r(f)ao 



(2.22) 



A plot of this quantity as a function of o"o is given in Figjsj It is clear that there is a minimum 
value of the distance between the two branes, the minimum occurs at ctq = and its value 
is 

_ 2^/2^r {\)^/K 
Amin — T~{^ (2.23) 

Since k is related to the world-volume gauge field, we see that only for zero electric field the 
two branes can annihilate. For large (Tq the distance A between the two branes grows linearly 



Figure 5: A for k = 1 in tlie Callan-Maldacena case 



with (jQ. We can now solve (2.22) for do by keeping fixed the distance between the branes A 



and the number of strings, which is done by keeping the charge parameter k fixed. We obtain 



-I 



A2 ± VA4 - 404^2 
2^^2 



(2.24) 
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2 /^p ( — \ 

where the numerical constant a is given by = — z-^gy^. There are two solutions which, for 
large A, behave as ^ 

QK A , . 

(70 ~ — , do ~ - 2.25 
A a 

In the first case, the "thin throat" branch, the radius of the throat goes to zero as A — t- oo. 
In the second case, the "thick throat" branch, the radius of the throat grows linearly with A 
as A — oo. 



3 Heating up DBI solutions by the blackfold approach 

Our aim in the rest of the paper is to turn on the temperature and analyze what happens 
to the Blon solution described in the previous section. We propose in this section that 
the appropriate framework in which to heat up DBI solutions is provided by the blackfold 
approach [17^ [T8] . We will motivate this here by first casting the DBI EOMs in a form that 



makes apparent the generalization to the thermally excited discussed in Section 13.1 

The latter involves the EM tensor for a black D3-F1 brane bound state, which is subsequently 
obtained in Section [3.21 



In Section 3.3 we give a detailed argument for why our thermal D-brane probe, based 
on the blackfold approach, should give a more accurate way to probe finite temperature 
backgrounds in string theory than the Euclidean DBI probe method that has been used in 
the literature. 

3.1 Extrinsic embedding equations from DBI action 

In this section we write a general expression of the EOMs for the DBI action. This takes 
the form of a set of extrinsic embedding equations. Below we shall use this to connect to the 
blackfold approach of [ITJ [18] and thence to go to the thermal case. 

Our starting poinl|^is the DBI action for the D3-brane given in (2.1). Before considering 



the EOMs we first obtain the world-volume EM tensor. This can be done by varying the 



action (2.1) with respect to the world- volume metric 7ab in (2.2), i.e. 



We find 

Td3 V-det(7 + 27r/2F) 



%/7 



((7 + 27rliF)-^r' + ((7 + 2<2i7)-i)f'a (3.2) 



where we defined the determinant 7 = — det(7a{,). 

We now consider the EOMs for the DBI action (2.1 ). These are found by variation of the 
embedding map X^((t"). We first notice that the Born-Infeld term in (2.1) only contributes 
through the variation of the world-volume metric. We compute 

= g^.u^xdaX^'dbX^dX^ + g^x{^aX^'^l,5X^ + d^X^dadX^) (3.3) 



Our considerations are easily generalized to Dp-branes 
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Hence we can write the variation of the Langrangian density of (2.1) as 
1 



+ 

giving the EOMs 

1 ■ — ab , 



4! 



3 nbcd 



e'''>"'daX''dbX''d,Xf{4dddX^C^,pX + ddX^C^,pa,xSX^) 



(3.4) 



VlT'''gf,,,xdaX>'dhX'' - db{^T-')g^xdaX^^ - ^T'^'g^x^udaX^d^X^ 
-^T-'g^^xdadbX^ + ^e-''''daX^d,X''d,XPdaX^{C^,p^,x 



We define now the projector h^'^ along the tangent directions to the D3-brane 

hf"" = -f^''^aX^'^bX'' 

along with the projector J-f^'^ along the orthogonal directions, defined as ±f^'^= g^'^ 
Using these we can define the extrinsic curvature tensor for the embedding 



(3.5) 

(3.6) 
- h^"'. 

(3.7) 



See [T7] for more on these geometrical quantities. We define furthermore the partial pullback 
of the RR five-form field strength F(5) = dC(4) 



-^Xabcd 

daX^^d.X^'d.XPddX'' 

and the D3-brane RR charge current 



abed 



\/7 



abed 



Projecting now the EOMs (|3^ with we can write the resulting EOMs as 



abed i 



Xabed 



(3.8) 



(3.9) 



(3.10) 



using the above definitions. This is the extrinsic equation for the D3-brane, with the EM 
tensor given by (3.2). As explained in |17] the equation (3.10) is basically the Newton's 
second law of brane mechanics, with T"* replacing the mass and Kab^ the acceleration of the 
point particle, while the right hand side generalizes the Lorentz force for a charged particlej^ 
The significance of rewriting the EOM's of the DBI action as (3.10) is that the generaliza- 
tion to a thermally excited D3-brane becomes more apparent. In fact, as we shall see below, 
in some sense all we need to do is to replace the EM tensor (3.2) with that of a thermally 
excited brane in (3.10) and then we have the EOMs for a thermally excited D3-brane. How- 
ever, in order to do this we have to shift to a regime in which the EM tensor of a thermally 
excited D3-brane is known: When we have a large number N of coincident D3-branes and 
gsN is largej^ In this regime the D3-branes backreact on the geometry of the background and 

^If we instead project the EOMs ( |3.5[ ) with h''^ we obtain the equation for conservation of the EM tensor 
Tab on the brane. See [T^ for general comments on tliis. 

^That QsN should be large follows from demanding that the curvature length scale of the supergravity 
solution of coincident extremal D3-branes, which is {NT-dzGY^'^ ~ {gaNY^'^lg, should be larger than the 
string length Is. 
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we can compute the EM tensor from the supergravity solutions of non-extremal black D3- 
branes. Doing this, we are essentially employing the recently developed blackfold approach 
for higher-dimensional black holes \17\ [TS] . 

In the blackfold approach \n\ [T5] one finds new black hole solutions by bending flat branes 
on curved embedding geometries. Indeed, as found in |17j . a brane with EM tensor Tab and 
current Jabcd with respect to the field strength F(5) would obey the extrinsic equation (3.10) 



to leading order in the regime where the brane thickness, as defined by its backreacion on 
the background, is small compared to the length scale of the curvature of the embedding 
submanifold defined by X^(a°'). This regime is parallel to the regime of the DBI action, i.e. 
the DBI action is the leading order action in the regime where the length scale of the open 
string excitations are small as compared to the length scale of the embedding geometry and 
of the variation of the gauge field on the brane. Higher order terms would for instance take 
into account derivatives of the Fat gauge field strength. 

Below we take the steps to write down the EOMs for thermally excited D3-branes. We first 
find the appropriate EM tensor in the next subsection, and remark also on the relation to the 
extremal case. Then we write down in Sec. [4] the correct EOMs and conservation equations in 
accordance with the blackfold approach [T71 [TH] and present the solution. We will also show 



how the extrinsic equations (3.10) for the thermally excited case can be obtained from an 



action, and discuss the relation of the latter to the DBI Hamiltonian. 

3.2 Energy-momentum tensor for black D3-F1 brane bound state 



In this section we consider the EM tensor of D3-branes. The EM tensor (3.2) corresponds to 
a single D3-brane at zero temperature with a gauge field strength Fab turned on. We work in 
the regime in which ^ab and Fab vary so slowly over the brane that their derivatives can be 
ignored. The aim of this section is to find the EM tensor of a D3-brane still with the gauge 
field strength turned on, but now at non-zero temperature and with a large number ^ 1 
of parallel stacked D3-branes and with ^^A^ large. 

More on DBI case 

Before turning to the non-zero temperature case we take another look at the zero-temperature 



EM tensor (3.2). Consider the special case in which the world- volume metric is diagonal, 



700 = — 1, and the only non-zero component of the two- form gauge field strength is Fqi. Then 



the EM tensor (3.2) becomes 



r°° = , = , = -7^1— Si= , = -y%3\/l -^2 i = 2,3 (3.11) 
where we defined 

E^^Foi (3.12) 
V7ii 

In the bulk, turning on the electric flux Fqi on the brane corresponds to having a number 
of F-strings along the direction of the brane. Call this number of F-strings k. From the 
world- volume point of view we can say that we have k units of electric flux. We can relate k 
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to Foi as 



. = -To3/ <<.w^^^H5pMS = ^/ ,,wv^-£= (3.13) 

If we consider the case in which jab and Fqi are constant the above equation can be written 
as kTpi = Vj_7b3-E'/Vl — E'^ where V± is the area in the a^'^ directions perpendicular to the 
F-strings. Using this we see that the mass density of the brane can be written as 



Too = ^T^, + ^ (3.14) 

We recognize this as the 1/2 BPS mass density formula for the D3-F1 brane bound state in 
the case of a single D3-brane and k F-strings. 

Energy-momentum tensor for D3-F1 bound state from black brane geometry 

We now turn to obtaining the EM tensor for D3-branes with an electric field on, corre- 
sponding to k units of electric flux, at non-zero temperature. We can obtain this from a black 
D3-F1 brane bound state geometry assuming that we are in the regime of large and gsN. 
The D3-F1 black brane bound state background has the string frame metric |21j 



ds'^ = D-^ H-^-fdt'^ + dxj) + H-^dxl + dxl) + D-^ H^f-^dr^ + r^dnl) (3.15) 
where 

/ = l-4, ^i^ ^osinh^« D-' = cos\ + sm^CH-' (3.16) 

and with dilaton field (j), Kalb-Ramond field -6(2)) and two- and four- form Ramond-Ramond 
gauge fields C(2) and C(4) given by 

g2</, = ^ = sinCiH-^ - 1) cotha 

(3.17) 

C23 = tan C{H-^D - 1) , C0123 = cos QD{H-^ - 1) coth a 

We now proceed to read off the EM tensor, the D3-brane and F-string currents, and the 
thermodynamical parameters as seen by an asymptotic observer. However, we first need to 
consider how the non-trivial world- volume metric jab can enter in this. Reading off the EM 



tensor from (3.15) we find it in the (t,x*) coordinates for which the world- volume metric is 



just -dt^ + Y,^{dx'Y . Instead we want the EM tensor with a world- volume metric of the form 

jabda^da^ = -dal + judal + 722^0-! + 733rfo-| (3.18) 

since this is the most general form needed for our computations in this paper, i.e. a diagonal 
metric without red-shift factor. To transform the resulting EM tensor in the (t, x*) coordinates 
to the above world- volume coordinates cr" we can simply make the rescaling t = cr'^, = 
^/Tiicr*, i = 1,2,3. One could infer here that such a rescaling is problematic since in general 
the world- volume metric jab varies according to where we are situated on the brane. However, 
to the order we are working in we are precisely suppressing the derivative of the metric in the 
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EM tensor, just as they are suppressed in the DBI EM tensor (3.2). With this in mind we 
read off the EM tensor in the a" world-volume coordinates 

= ^r^3ro^(5 + 4sinh2a) , T^^ = - "^T^.rl {1 + 4 sinh' a) 



(3.19) 



using for example [22]. The D3-brane current is 



27r^r^ 

j0123 ^ D3 ^^4 ^^gj^ ^ g-^j^ ^ (^3 20) 



The number of D3-branes in the bound state is A^. Thus, using (3.20) we find 



'D3 

Furthermore, imposing that we have k F-strings gives 



N 

cos C^'o '^osh a sinh Q = — 5— — (3-21) 



N ~ Tf 



- [ da^da^^/^^^tanC (3.22) 

1 JV'2:i 



'V23 

We also give the thermodynamic quantities associated to the horizon, that will be used in 
what follows. The temperature T and entropy density 5 are 

T = ^— — , 5 = 27^3^^oV.o^r^cosha (3.23) 

vrrocosha U6 (d) \ ; 

while the local D3-brane and F-string chemical potentials are 

//[^g'^'^''' = tanh a cos C , /ip°'^'^^^ = tanh a sin C (3.24) 

Extremal limit 

As a check we consider here the extremal limit of the EM tensor of the black D3-F1 brane 



bound state. We see from (3.21) that the extremal limit is to take tq — )• keeping C and 
Tq cosh Q sinh a fixed. This gives T = as it should and the EM tensor 

^00 ^ NTm ^ 5.11 ^ _^ii ^ ^ _Y^NTj,s cos C, i = 2, 3 (3.25) 

cos C COS C 



along with the formula (3.22) for k. We see that the above formulas match (3.11) and (3.13) 



provided we identify E = sm( and put A*" = 1. 

3.3 Argument for new approach to thermal D-brane probes 

As reviewed in the introduction, a number of papers in the literature (see |23t [131 HM [T5] and 
later works) have used the classical DBI action to probe finite temperature backgrounds in 
string theory. In short, this method consists in Wick rotating both the background as well 
as the classical DBI action, then finding solutions of the EOMs from the classical Euclidean 
DBI action and finally identifying the radius of the thermal circle of the background with the 
radius of the thermal circle in the Euclidean DBI action. From the classical solution one can 
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then evaluate physical quantities for the probe such as the energy, entropy and free energy. 
We dub here this method the "Euclidean DBI probe" method. 

We give here a detailed argument for why our thermal D-brane probe, based on the 
blackfold approach, should give a more accurate way to probe finite temperature backgrounds 
in string theory than the Euclidean DBI probe methodj^ 

We begin by considering in more detail the Euclidean DBI probe method. For simplicity 
we stick to a D3-brane (thus in type IIB string theory) but our considerations apply to any 
Dp-brane. Consider a type IIB string theory background with metric g^u and with a RR five- 
form field strength ^(5) turned on (for simplicity we do not consider other RR field strengths 
and we also assume a constant dilaton). As shown above, the EOMs for the D3-brane DBI 
action (2.1) take the form 

T^BlKab' =^'^ ^r'^'^Fxabcd (3-26) 



where Tggj is given by ([O) and J"^^'^ by S. We now perform a Wick rotation on the 
background t = He where t is the time coordinate of the background and Ie the corresponding 
direction in the Euclidean section of the background. Similarly we also perform a Wick 
rotation for the world- volume time r = ite- Then the EOMs (3.26) become 

{TErB^l{KE)ab' = {^Er^liJEr'"'{FE)xabcd (3.27) 

where the subscript E means that it is the Wick rotated quantity where the Wick rotation in 
both the bulk and on the world-volume are treated as simple linear transformations on the 
tensors, e.g. {Te)i^^i = — ^dbi ^^'^ s° ^'^^ easily shown that one also obtains the 

equations ( |3.27 ) as the EOMs of the Euclidean DBI probe in the Wick rotated background, i.e. 
by varying the Wick rotated DBI action in the Wick rotated background. We can conclude 
from this that there is a one-to-one map between solutions of the EOMs for a DBI probe 
in a thermal background and the solutions of the EOMs for a Euclidean DBI probe in the 
Wick-rotated thermal background. 

To solve EOMs corresponds to solving certain differential equations under the restriction 
of certain boundary conditions. The above equivalence between solving the EOMs for a DBI 
probe in a thermal background and the EOMs for an Euclidean DBI probe in the Wick- 
rotated thermal background only means that the differential equations are the same, instead 
the boundary conditions are different. Thus, the equivalence works only locally. Instead there 
are global differences in being in the Wick-rotated frame or not since one imposes different 
boundary conditions, in particular regarding the thermal circle direction tE. In the Euclidean 
probe method one uses the Wick rotated version of the usual static gauge Ie = te and 
one imposes that the size of the thermal circle is the same for the probe as for the (Wick 
rotated) thermal background. Thanks to the static gauge this boundary condition does not 
enter in the EOMs for the Euclidean probe and is in that sense a global condition on the 
solution. Therefore, one could think that this global condition is enough to ensure that the 
probe, using the Euclidean DBI probe method, is in thermal equilibrium with the background. 
However, we shall now argue that this global condition is not enough since the requirement of 
thermal equilibrium between the probe and the background changes the EOMs of the probe 



''The considerations of this section were independently worked out by Roberto Emparan. 
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by changing the EM tensor, and hence the requirement of thermal equihbrium changes the 
probe not only globally but also locally. 

Consider the example of a D3-brane with zero world-volume field strength Fab = in the 
AdSs X background. In the probe approximation this is described by the DBI action. The 
EOMs are of the form ( |3.26[ ) with Tg^j = -Tdsj"^. We notice that the EM tensor locally is 
Lorentz invariant. This conforms with the fact that the electromagnetic field on the D3-brane 
is in the vacuum state. 

We now turn on the temperature in the background. Thus, the background is either hot 
AdS space, or an AdS black hole, depending on the temperature, times the S^. The D3- 
brane should be in thermal equilibrium with the background. This means in particular that 
the DOFs living on the D3-brane get thermally excited, acquiring the temperature of the 
background. Among the DOFs are the ones described by the electromagnetic field Fab living 
on the brane. The quantum excitations of Fab are described by the maximally supersymmetric 
quantum electrodynamics locally on the brane (since one can see by expanding the DBI action 
that one locally has Maxwell electrodynamics for small Fab which then is supplemented by 
the superpartners from the Af = 4: supersymmetry) . This means that for small temperatures, 
near extremality, the EM tensor takes the form of the Lorentz invariant piece plus the EM 
tensor corresponding to a gas of photons and their superpartners. Consider a particular point 
q on the brane. We can always transform the coordinates locally so that jab = Vab at that 
point. Then the EM tensor at q takes the form 

Tab = -TmVab + riJ^'^^ , T^T'^ = P , T^r^ = p , i = 1, 2, 3 (3.28) 

where T^^^ is the contribution due to the gas of photons and superpartners, having the 
equation of state p = 3p = 7r^r^/2 (the power of follows from the fact that = 4 
supersymmetric quantum electrodynamics is conformally invariant). Thus, we see that the 
fact that we have local DOFs living on the D3-brane means that the EM tensor is changed once 
we turn on the temperature. In terms of the EOMs for the probe we see that they are given 
by (3.10) with the EM tensor ( 3.28| ) (one can easily find this EM tensor for general world- 
volume metric Jab)- Therefore, the EOMs are clearly not the same as those of (3.26) where 
-^DBi ~ ~Tb37"'^- Indeed, the EM tensor (3.28) is no longer locally Lorentz invariant, which 
is in agreement with the fact that the brane has an excited gas of photons and superpartners 
on it. 

In conclusion, the above example clearly illustrates that the requirement of thermal equi- 
librium affects the probe not only globally but also locally in that the EOMs change from 
those given from the DBI action. This shows that the thermal D-brane probe is not accurately 
described by the Euclidean DBI probe method}^ 

Finally, we note that our description of the thermal D-brane probe is in accordance with 
the above example. Indeed, if we expand the EM tensor for a non-extremal D3-brane 

2 2 

r°° = Yr23r^(5 + 4sinh2a) , = -y*yr23r^(l + 4sinh2 a) , i = 1,2,3 (3.29) 



Note that the difference between the thermal D-brane probe and the Euchdean DBI probe is not due to 
backreaction. A backreaction would mean that Kab'' should change. Instead, demanding thermal equilibrium 
means that the EM tensor changes. Thus, even in the probe approximation the Euclidean DBI probe is not 
accurate. 
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for small temperatures we get 

37r2 



NTr,3 + — iV^r^ , = f'{-NTD3 + —N'^T'') , i = 1,2,3 



vr 



2rp4\ 



(3.30) 



We see that this is precisely of the form (3.28) (with an extra factor of N in the leading part) 
with the leading part being locally Lorentz invariant and the correction corresponding to a 
gas of gluons and their superpartners (with the usual factor of 3/4 since we are at strong 
coupling pi]). This is thus in accordance with our general arguments above. Finally, we 
note that whereas in the regime of validity of the DBI action the near-extremal correction 



to the EM tensor (3.28) should be computed quantum mechanically, we can use classical 



supergravity to compute the full EM tensor at finite temperature in the large N and large 
gsN regime since the classical approximation is reliable in this regime. 



4 Thermal D3-brane configuration with electric fiux ending in 
throat 

In this section we study the D3-F1 configuration at finite temperature in hot fiat space. 
We derive the EOMs for the embedding directly and supplement this with an action-based 
derivation. We then proceed by solving the equations and comment on the relation to the 
DBI solution and Hamiltonian. Finally, we discuss the regime of validity of the solution. 



4.1 D3-F1 extrinsic blackfold equation 

Our setup is specified by exactly the same type of embedding and boundary conditions as 



discussed in Sec. 2.1 for the extremal case, see Eqs. (2.3)-(2.5), (2.7), (2.8) and the illustration 



in Fig. [T] With the results of the previous section, we are thus ready to compute the D3-F1 



extrinsic blackfold equation (3.10), where the right-hand side is zero since our background is 



lOD Minkowski space-time so there is no five-form field strength. For the left-hand side we 



need to compute the extrinsic curvature tensor (3.7) for the embedding described in (2.3)-(2.5). 



The resulting non- vanishing components are given by 



a sm 



l + z'(a)2 ' 



'■22 



l + z'(a) 



'■33 



l + z'{ay 



(4.1) 



with i = 1,2,3. Since Kl^ and K^^^ are proportional, the EOMs (3.10) for this case becomes 
simply T'^'^Kab'^ = 0. Here Tab is the EM tensor (13.19^ of the black D3-F1 brane system, 



where the collective coordinates ro, a, C, are now promoted to be functions of a. We then find 
the following EOM for the D3-F1 blackfold 



z'{l + z' 



2 1 + 4 cos'' C sinh'' a 
a 1 + 4 sinh^ a 



(4.2) 



The precise form of the functions C(c) and a{a) entering this equation (as well as ^0(0")) follow 
from a number of constraints on the solution, as we will now discuss. 
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Constraints on solution 



Beyond the EOM (4.2) we also have further constraints that follow from charge conservation, 



of both D3-brane and F-string charge, and constancy of the temperature. From the point 
of view of the general blackfold construction discussed in Ref. [E] the constancy of the 
temperature and angular velocity is a consequence of stationarityj^ Since in our case the 
blackfold is static only the temperature is relevant. Furthermore, charge conservation can 
be seen to follow from additional EOMs, as will be explained and derived in detail in the 
forthcoming paper Ref. [19| which generalizes the construction of neutral blackfolds to charged 
blackfolds. 



We start by considering the F-string charge, given in eq. (3.22 ), which should be conserved 



along the a direction. Inserting the induced metric (2.5), this gives that 



where the constant k is defined af 



10 



a tan C 



in terms of the conserved charges and k. Note that we can also write (4.3) as 

cos C - 



1 



(4.3) 



(4.4) 



(4.5) 



1 + S 



We also need to ensure the conservation of F-string chemical potential |19j (the quantity 
defined as / da 



Tii/^Fi'^'^^^ with ^p"^^^^ the local F-string chemical potential given in (3.24)) 



in the world- volume directions transverse to a. This is automatic since we impose spherical 
symmetry for the two-spheres parameterized by the coordinates 0, (/> transverse to the F-string 
direction a on the world-volume of the D3-brane. 



The equation (4.5) ensures F-string charge conservation by conserving the number of F- 



strings k assuming the conservation of the number of D3-branes A^. The latter is imposed 



using the formula given in eq. (3.21), while constancy of the temperature is imposed using 



Eq. (3.23). Thus, from Eqs. (3.21), (3.23) and (4.5) we can eliminate ro(o") and ("(cr) yielding 

1 + ^ . (4.6) 



the constraint 



sinh a vr^ NT^ 



cosh'^ a 2 Td3 V a"" 
Below we solve this equation explicitly for cosh a. 

Summarizing, we note that for given (T, A^, K) the cr-dependence of the two functions C(c) 



and a((7) is determined by the two equations (4.5) and (4.6) and. Given the solution of a(a"). 



one can then compute the thickness ro((T) from (3.23) 



^More precisely, in the blackfold approach the intrinsic EOMs are conservation of the world volume EM 
tensor. For stationary solutions to these one finds that the temperature and angular velocities (if present) are 



constant 

10 



Note that, in a slight (but meaningful) abuse of notation this differs by a factor of 1/A'' from k used in the 
DBI analysis in Sec. [2j 
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Action derivation 



It is also instructive and useful to consider an action derivation of the EOMs (4.2). In fact, 



as we shall see later, the action derived below helps in finding an analytic solution of the 
equations of motion. Moreover, it can be viewed as the action that replaces the DBI action 
when thermally exciting the D3-F1 system in the regime of large number of D3-branes and 
large QsN. 

To write down the action, we use the fact that Refs. \n\ IT^ showed that for stationary 
blackfolds the extrinsic blackfold equations can be integrated to an action which is proportional 
to the Gibbs free energyp] We therefore use the thermodynamic actior 



F=M-TS 



(4.7) 



where J- = J-{T, N, k) is the free energy appropriate for the ensemble where the temperature 
T and number of D3-branes A'^ and F-strings k are fixed. Here the total mass M and entropy 
S are found by integrating the energy density T^^ in (3.19) and the entropy density S in 



(3.23) over the D3-brane worldvolume. 

M = Y^DS I dV^3)r^{5 + 4 sinh^ a) , S = l^^T^^ j dV^^,)rl cosh a 



(4.8) 



These are thus regarded here as functionals of the embedding function -z(ct), that give the 



actual total mass and entropy of the system when evaluated on-shell. Using (4.^ 
then find the action functional 

^2 



in (4.7) we 



IT' 

~2 



T^3 J dV(3)r^(l + 4sinh2 



a 



Eliminating ro{a) using (3.23) and using the induced metric (2.5) we get 

/ da^/l + z'{cjYF{a) 

J an 



T 



2T2 



where we introduced the function 



a 



> 4 cosh^ a — 3 
cosh^ a 



(4.9) 



(4.10) 



(4.11) 



Note that we integrate from ctq to infinity in (4.10) according to the boundary conditions 
discussed above. 



The function F{a) defined in (4.11) is a specific function of a for given T, N and k, as 



seen from (4.6). This means that when we vary the action (4.10) with respect to z{a) the 



function F{a) does not vary. Performing this variation in the action (4.10) we then find that 
the EOMs take the iovu^ 

( z'{a)F{a) V , 



(4.12) 



^^This was done in Ref. [T^ for neutral stationary blackfolds and generalized to the charged case in Ref. |19] . 
Note also that this implies that for stationary blackfolds the extrinsic equations are equivalent to requiring the 
first law of thermodynamics. 

^^More properly, the action is / = PF, but since /3 = l/T is constant we directly use the free energy T. 

^''This is computed as + z''^F) = ^ 5z' and adding a total derivative to the last formula. 
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which we wih use below when solving the system. It is noteworthy that the EOM (4.12) is 



exactly as that of (2.14) for the DBI case with -Fdbi replaced by F in (4.11). The same holds 



in fact when comparing the free energy (4.10) and the DBI Hamiltonian (2.13). In Sec. 4.3 
we will further comment on their relation. 



As a check on the action approach, we now show that (4.12) obtained from the thermo- 



dynamic action is consistent with the D3-F1 extrinsic blackfold equation (4.2). From (4.12) 
we find 



z'{l + z 



I2\ 



F^ 



(4.13) 



Computing F'{a) from (4.11) and using ( |4.5[ ) and (4.6) to eliminate a'{a) and k, we get 

2a 



F' 



cosh'^ a 



1 + 4 cos^ C sinh^ a 



(4.14) 



Using this with (4.11) we indeed obtain (4.2) from (4.12). Note that the equivalence of these 



two equations implies that the blackfold EOMs are equivalent to requiring the first law of 
thermodynamics, since the latter follows from extremizing the free energy functional. 

4.2 Solution and bounds 



We now proceed solving the EOM (4.12) subject to the constraints (4.5) and (4.6). The latter 



imply a bound on the temperature as well as the world-volume coordinate a parameterizing 
the size of the two-sphere, which we will first discuss. Then we turn to the explicit solution of 



the constraint (4.6) in terms of the constants (T, N, k) and subsequently present the solution 



of the EOMs for the profile z{(t) 



Bounds on the temperature and a 

Considering the left hand side of the constraint ( |4.6[ ) it is easy to check that sinh a / cosh'^ a 
is bounded from above, with a maximal value 2\/3/9 for cosh^ a = 3/2. Using (4.5) we thus 
get an upper bound for the temperature 



< Tbnd cos C 



(4.15) 



for given (A^, A;) and a, where we defined the temperature 

/4V3Td3\' 



nd 



I 97r2Ar 



(4.16) 



From (4.15) it is obvious that we furthermore have the weaker upper bound T < Tbnd which 



only depends on N and thus not on k and a. The temperature Tbnd is the maximal temper- 
ature for N coincident non-extremal D3-branes, which can never be reached in the presence 
of non-zero F-string charge. For future convenience we furthermore define the rescaled tem- 
perature 

f^-^ (4.17) 

^hnd 



which will simplify expressions below. Note from (4.15) that we have T < 1. 
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Since a takes values in [(To,oo) we see that cosC in (4.5) is minimized for a = ao- It then 
follows that the upperbound ( 4.15| ) can be stated more accurately as 

T < Tbnd ( 1 + — ) (4.18) 



cr?, 



for given N, k and ao of the brane profile z{a). Note however that generically this bound can- 
not be saturated for a given N, k and do, so that generically the maximal temperature is lower 



than this bound. Finally, we remark that using (4.5) the upper bound on the temperature 



(4.15) can be turned into a lower bound on a 
for given (T, k). 



k(T- 



1) 



(4.19) 



Solving cosh a in terms of a 



We now present the explicit solution of the constraint HM. Using the upper bound (4.15) 



on the temperature along with (4.5) we can define the angle S{a) by 



(4.20) 



where we restrict the angle to be in the interval < 5{a) < tt /2. In terms of this angle the 



constraint (4.6) can be written as 

4 cos^ 5 



27 



cosh a — cosh a + 1 = 



(4.21) 



which is a cubic equation in cosh^ a. Thus, it has three independent solutions for cosh^ a in 
terms of 5. The first solution is 



cosh^ a 



3 cos I + \/3 sin 



(4.22) 



2 cos 5 

As 5 goes from to '7r/2, the right hand side increases monotonically from 3/2 to infinity. 
Substituting (5 — )■ — (5 in (4.22) we find the second solution 

, 9 3 cos I — a/S sin I ^ ^ 

cosh^ a = ^ — ^ (4.23) 



2 cos 5 

Here the right hand side decreases monotonically from 3/2 to 1 as (5 goes from to 7r/2. 
Finally, by substituting 5 



(5 — 27r in (|4.22|) we find the third solution 
cosh^ a 



3cos^ + \/3sin^ 



(4.24) 



2 cos 6 

This solution can be immediately discarded since it decreases from — 3 to — cxd as 5 goes from 
to 7r/2, and a has to be real. 



Turning to the two solutions (4.22) and (4.23) we note that both of them respect that 



cosh^ a > 1. In the extremal limit one takes a — )• oo. Therefore, the solution branch connected 
to the extremal solution is the first solution (4.22). For this branch the energy density {T^^ in 



(3.19)) for each value of o" increases as the temperature T increases holding (A^, k) fixed. Thus, 



this is the thermodynamically stable branch with positive heat capacity. Instead in the second 



solution branch (4.23) the energy density at each a decreases with increasing temperature. 



thus resulting in a negative heat capacity. This branch is connected to the neutral 3-brane 
with a = 0. The two branches meet at the point cosh^ a = 3/2. 
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The solution 



Finally, we turn to the solution of the EOM (4.12). Imposing the boundary condition (2.8) 
we find that the solution takes fornix 



Imposing the boundary condition (2.7) we find 
z{a) = I da' { - 1 



F{a) 



a 



, 4 cosh^ a — 3 
cosh^ a 



(4.25) 



(4.26) 



with a > ao and where we repeated the definition of F{a) from (4.11) given in terms of 
cosh^ a. Furthermore, cosh'^ 0(0") is given by one of the two branches (4.22) and (4.23). The 



result (4.26) describing the heated up Blon solution is one of the central results of this paper. 



Just as the Blon solution, it is obtained in a probe approximation, which is discussed in more 
detail in Sec. l44l 

It should be emphasized that for each of the two branches we have an explicit expression 
for the derivative of the brane profile z'{a) so that the mass and entropy for each solution 



branch can be obtained explicitly from the integrands in (4.8), (4.10). This gives 

F(a) 



M 



2T2 



vr 



rp4 



da — , 

-0 VfW) 



s 



2T2 
2^5 



vr 



da — ^ 



F^iao) 
F{a) 



2 4 cosh a + 1 
cosh^ a 



: a 



F2(cro) cosh^a 



We also give the integrated chemical potentials [19j that follow from (3.24 
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/UD3 = 47rrD3 



da 



TO 



da- 



^F^{a)-F^ao) 
F(a) 



cr^ tanh a cos C 



tanh a sin C 



,.0 ^F\a)-F^iaoj 
These satisfy the first law of thermodynamics and Smarr relation 

dM = TdS + HDsdN + ^ipidk , 4(M - /xds^V - fiFik) = 5TS 



(4.27) 
(4.28) 

(4.29) 
(4.30) 

(4.31) 



The integrations in (4.27)-(4.30) can be performed numerically or analytically in certain limits. 



The solution presented above gives the profile of a configuration of N coincident infinitely 
extended D3-branes with k units of F-string charge, ending in a throat with minimal radius 
ao, at temperature T. The configuration is sketched in Fig. [2j As explained in Sec. 2.1 



we can 



construct a corresponding wormhole solution by attaching a mirror solution as illustrated in 
Fig. |4j This will be discussed in Sec. [5j In the next subsection we examine two limits of the 
branch ( |4.22 ) connected to the extremal configuration, which will be used in Sec. [5] App. |A] 



presents some features of the branch (4.23) connected to the neutral 3-brane, which will not 
be further discussed in the main text. 



^■'Notice that it follows from the expression (4.251 that the derivatives d"a/dz" of the inverse function o-{z) 
vanish for odd n at z = z(o"o). This is a necessary requisite for the smoothness of the wormhole solution of 
Section [5] 

^^We redefine these for convenience by an extra factor of the D3-brane and F-string tension respectively. 



22 



4.3 Analysis of branch connected to extremal configuration 



Our main focus in this paper is the solution (4.26) for the branch (4.22) connected to the 



extremal configuration, which we start analyzing here by considering two physically relevant 
limits. We then go on to showing the equivalence of our thermal action (free energy) ( |4.10 ) 
in the zero temperature limit to the DBI Hamiltonian. 

One interesting limit is the limit of small temperature, which enables to see in a small 
temperature expansion the effect of heating up the Blon as compared to the extremal case. 



The temperature enters the solution through the function a{a), which is given by (4.22) for 



the branch connected to the extremal configuration. From (4.20) it is seen that for small 
temperatures 6 is close to tt/2, so that 



cosh a 



3^3 
2 cos 6 



V3 



cos (5 cos^ 6 + 0{cos^ 5) 



12 27 

This gives for the function F{a) defined in ( 4.11| ) the expansion 



COS 6 



cos^ 6 



cos3 5 + 0(cos^ 5) 



(4.32) 



(4.33) 



9 27 243 

Keeping the first two terms in the expansion we then find from the solution (4.25) the result 



-z'{a) 
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1 + 



'1 + 



a" 



(4.34) 



Taking the zero temperature limit, we see that the first term agrees with the Blon solution 



(2.15). The second term describes the leading order effect due to heating up the Blon. Note 
that the factor in front of does not blow up as a — )• (Tq. This expression is a good 
approximation for all a > ao when (T, A^, k) and (Tq are given such that T^y^l + k'^/cq <^ 1. 

Another interesting limit is a j ^fk ^ 1. This limit describes the profile for large cr, i.e. 
in the region close to the fiat D3-brane at infinite o. In this limit we have to leading order 
cos 5 = T^. Therefore, from ([O2I) we see cosh a is a function of T only. For large a/y/H we 
thus have 

F((t) = (T'^g{f) + ©(k/cj^) (4.35) 
where g{T) is a function that increases from to 4/3 as T goes from to 1. For sufficiently 



large a we find therefore from (4.26) the behavior 



9{T) (y 



(4.36) 



If we require (t/^/k ^ 1 for all o", the limit corresponds to aQ/^/H » 1, which can be regarded 
as the limit in which the F-string charge is taken to be small. In this case it follows from 
(|4.35D that F(cr)/F(ao) ~ a'^/al and hence 



0"n 



(4.37) 



which agrees with (2.15) in the limit gq/^/k ^ 1. 
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Comparison with DBI Hamiltonian 

Finally, we point out the relation of the DBI action with our thermal D3-F1 brane action (free 



energy) in (4.10). To this end we compute the T — >• limit of our action. In particular, using 



the small temperature expansion (4.33) for the branch connected to the extremal configuration 



along with the expression (4.20) for cos 5 and the definition (4.16) we find by inserting in the 



action (4.10) that 



lim T 

T->0 



DBl|To3^iVTo3 



(4.38) 



This shows that our thermodynamic action (4.7) can be viewed as the thermalization of the 
DBI Hamiltonian (2.13). The replacement Td3 — t- NTq^ induces an extra factor of N in front 



of the DBI Hamiltonian (2.13), and at the same time includes an extra factor of 1/A^ in k, 



yielding the one defined in (4.4). It is satisfying to recover these multiplicative factors of 



N since the DBI action is valid for a single D3-brane only and our thermal system contains 
instead N D3-branes. Moreover, this relation implies that the systems are properly connected 



also off-shell and makes manifest that the extremal limit of the branch (4.22) of our solution 
gives the Blon solution reviewed in Section [2] 

4.4 Validity of the probe approximation 



In the blackfold approach used to obtain the EOM (4.2), the D3-F1 system is considered in 



the probe approximation. We therefore need to determine the conditions for which this is a 
valid approximation. If we denote the transverse size scale of the D3-F1 geometry by rg, there 
are two conditions that have to be met across the embedded surface in order for gravitational 
backreaction to be negligible: 

• The size scale should be much less than the 5^ radius a of the induced geometry: rg <C cr 



The size scale should be much less than the length scale 
of the induced geometry: rg <C -^^curv 



of the extrinsic curvature 



The transverse size scale rg of the D3-F1 geometry (3.15) is determined by the harmonic 



functions /, H and D in (3.16). The three corresponding scales are ro, tq sinh^/^ a and 



ro(sinCsinha)^/^ respectively. Since sin^ < 1, we only need to consider the first two scales 



and we let Vg be the largest of the two. For the main solution branch (4.22) connected to 
the extremal configuration the largest scale is then Tg = Tq sinh'^ a, since 1/2 < sinh^ q < oo. 
Moreover, it is in fact sufficient to consider the charge radius = Tq sinh a cosh a instead. 
Near extremality this is obviously true, but, more generally we have on this branch 1/3 < 
tanh^ a < 1, so that differs by a factor of order 1 from rg. For the other branch ( |4.23 ) 
connected to the neutral 3-brane, where we have < sinh^ a < 1/2, similar reasoning implies 
that it is sufficient to consider the radius = vq. Focusing on the main branch, the above 
conditions can thus be written as 



rc(cj) < CJ , rc(fT) < Lcurv(o-) 



where, using (3.21) and (4.5), the charge radius is found to be 



4 



Ty)3 



1 + 



(4.39) 



(4.40) 
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Considering the first condition rc{cr) ^ cr in (4.39) we see that rc{cr) is largest for a = ctq thus 
to satisfy rc(o") ^ a for all cr > fio it is sufficient to demand rc{(TQ) <^ ctq. 

To understand the second condition rc((7) ^ u in ( 4.39| ) we first need to determine 
Lcurv(o"). For this we need to compute the mean curvature vector Rf = 'y°'^K^f^. The curva- 
ture size is then given by Lcurv = with K 

of the embedding surface of the three-brane. In the case at hand we have 



K^Up, where Up is the unit normal vector 



1 



:(0,-z'(a), 0,0, 1,05) 



(4.41) 



in the coordinates {t, r, 6, (p, Xi) used in (2.3). To see this note that the tangent vector to the 3- 
brane is V = , ^ ^(0, 1, 0, 0, z'ia), O5), so that UpV = 0. Using the second fundamental 

tensor computed in (4.1) we then find after some algebra 



K 



1 



(1 + z'2)3/2 



z" + 2-{l + z 
a 



/2\ 



_ Fo 


'F' 


2" 


F 


~F 


a 



2 Fo 4 sinh^ a 



sin^ C 



(4.42) 



fT F 1 + 4 sinh^ a 

where Fq = F{ao)- To examine this further, we note that for the main branch we have 
sinh^ ce > 1/2 so that the the third factor in (4.42) is of order one. We have furthermore 
checked that the extrinsic curvature (4.42) reaches its maximum Xmax in the region a 



close to the end of the throat. We thus get 

K 



sin^ C 



(4.43) 



Hence the second condition in (4.39) requires to be in the regime where rc(o"o) ^ ao/sin^ C- 
However, this is already ensured by the stronger condition rc((To) <^ ao coming from the first 



condition in (4.39). Thus, in conclusion, it is sufficient to impose the condition 

rc(ao) < ao (4.44) 
to ensure the validity of the probe approximation. To get a better understanding of the 



condition (4.44), we now consider it in the regime where (Tq/vk is very small. Using the 
definition of k in (4.4 ) the condition becomes ^ Vkggls- It is interesting to note that the A^- 



dependence has canceled out in this condition. It is intuitively clear that the larger the number 
of F-strings, the greater the minimum radius should be in order to neglect backreaction. 
Moreover, for sufficiently weak string coupling it is always possible to be in the correct range. 



5 Separation between branes and anti-branes in wormhole so- 
lution 

5.1 Brane-antibrane wormhole solution 

In Section [4] we found a configuration where N coincident thermally excited D3-branes have 
an electric flux such that the equipotential surfaces of the electric field are on two-spheres. 



The shape of the configuration is given by the function z{a) in Eq. (4.26). As a ranges 
from 00 to (To the bulk coordinate z{(t) ranges from zero to z{ao) where the tangent of the 
brane is orthogonal to the tangent of the brane at cr = 00. Thus, we have k electric flux 
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lines all pointing towards a center, and as a decreases the density of the electric flux on the 
equipotential two-spheres increases. However, rather than having a central singularity as in 
the linear theory of Maxwell electrodynamics, the non-linear nature of the DBI theory on the 
D3-brane gives instead a bending of the brane preventing us from reaching the point where 
the singularity should have been. The spherically symmetric electric flux thus causes the 
bending of the D3-brane in the bulk creating a throat on the brane ending at r = ao and 
z = z{ao) = A/2 as illustrated in Figs. [l] and [2| As stated in Section|4]we should distinguish 
between solutions z{ao) = oo and z{ao) < oo. For z{ao) = oo we have an infinite spike. 
Instead, for z{aQ) < oo the solution as it is violates the conservation of D3-brane and F- 
string charge. To remedy this one can attach a mirror of the solution, reflected in xi around 
xi = z{ao) = A/2, as illustrated in Fig. |4j This we call a brane- antibrane wormhole solution. 

In this section we shall see that it is not possible to find finite-temperature solutions with 
z{ao) = oo, i.e. the infinite spike, unlike in the zero temperature case (see Section]!]). For this 
reason, we focus here on the brane-antibrane wormhole solution. However, a different way 
of generalizing the extremal infinite spike solution to finite temperature will be considered in 



In more detail, the brane-antibrane wormhole solution has N coincident D3-branes ex- 
tending to infinity at xi = 0, and anti D3-branes {i.e. oppositely charged) extending to 
infinity at xi = A = 2z((To). Then centered around r = we have a wormhole where for 
each constant xi slice you have a two-sphere where a flux of k F-strings is going through in 
the positive xi direction. At xi = z{ao) the two-sphere has the minimal radius r = ctq. The 
solution has altogether four parameters: The D3-brane charge N, the F-string charge k, the 
temperature T and the minimal radius ctq. It is important to note here that the extensive 



quantities given in (4.27)-(4.30) should be multiplied by a factor of 2 since we add a mirror 
of the solution. 

In this section we consider the separation distance A = 2z{ao) between the A^ D3-branes 
and A^ anti D3-branes for a given brane-antibrane wormhole configuration defined by the four 



parameters A^, k, T and ao. We consider only the branch (4.22) connected to the extremal 



solution. Thus, we have from (4.26) 



A 



(5.1) 



with F{a) given by (4.11), (4.22) and (4.20). From considering (4.11), (4.22) and (4.20) we 



see that F{a) depends on A^, k and T only through the variables k and T defined in (4.4) and 
(4.17). Therefore, the separation distance is a function with dependence A = A{T,ao, k). 



However, changing n corresponds to a uniform scaling of the system. Indeed, it is easy to 
show that A{T, (Tq, k) has the scaling property 



(5.2) 



Thus, it is enough to flnd A for k = 1. In detail, given A for a certain ctq and T with k = 1 
the general k conflguration is found by rescaling A — t- y^A and ctq — s- \/k(To while keeping T 
fixed. 
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5.2 Diagrams for separation distance A versus minimal radius (Jq 

In the following we examine the behavior of the separation distance A(T, o"o,k = 1) = 
A{T,ao). We set k = 1 in the rest of this section since, as argued above, one can always 
reinstate k by a trivial rescaling. 

We shall examine the separation distance A{T, ao) by using both numerical and analytical 
methods. We first evaluate A{T,ao) numerically. The behavior of A versus ctq for four 
different values of T is shown in Figures [g] and [Tj where we compare it with the case of zero 
temperature (Figjs]). We have chosen the values f = 0.05, 0.4, 0.7, 0.8. Note that < T < 1. 



A A 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

Figure 6: On the figures the solid red line is A versus (Jq either for T = 0.05 (left figure) or T = 0.4 
(right figure) while the blue dashed line corresponds to T = 0. We have set k = 1. 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Figure 7: On the figures the solid red line is A versus (Jq either for T — 0.7 (left figure) or T = 0.8 
(right figure) while the blue dashed line corresponds to T = 0. We have set /t = 1. 

By looking at the curves in Figures [6] and [7] we can see that there are new interesting 
features with respect to the zero temperature case displayed in Figure [S] We first notice that 
(To is bounded from below. Indeed, we found in Section |4] the bound ( |4.I9| ) 

(JO > amin = _ ^8)1/4 (^-3) 

This is in contrast to the zero temperature case where one can take ctq — >• corresponding 
to the infinite spike solution. For most values of T we further have the feature that when 
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increasing ao from <Tmin then A increases until it reaches a local maximum denoted by Amax- 
Increasing do further A decreases until it reaches a local minimum that we denote Amin 
(for most values of T this is also the global minimum). Increasing ctq further A increases 
monotonically and follows increasingly closely the zero temperature value of A as a function 
of ctq. In Figure |8] we have displayed the behavior of A at ctq = o"miin \na.x and Amm for all 
values of T. We see from this that there is a critical value of T given by T;, ~ 0.8 beyond 
which Ajnax and Amin ceases to exist. Note also that just before T reaches Tf, one has that 
A at o"o = cJinin is smaller than Amim unlike for lower values of T where Amin is the global 
minimum. 



A 




Figure 8: Amax (blue curve), Amin (red curve) , A at (J min (black curve) as a function of the 
temperature T for k = 1. 

We can conclude from the above that the there is no direct analogue of the infinite spike 



solution with ctq = and z{aQ) = oo for non-zero temperature. As one can see from (5.3) fio 
can not reach zero for non-zero temperature. Moreover, for T <Th the local maximum Amax 
is always finite, as one can see from Figure [8| Thus, there are no finite values of ctq for which 
A is infinite. However, in our paper [2Qj we examine a different way to find an analogue of 
the infinite spike solution at finite temperature. 

The plots in Figs. [6j [7] and [8] illustrate a number of interesting features: 

• We observe that for any temperature, for sufficiently large ao, the finite temperature 
curve is increasingly close to the corresponding extremal curve as ctq increases. In 
particular, in accordance with expectations, as the temperature is lowered a larger part 
of each of the two curves is close to each other. However, for any non-zero temperature 
T < Tb there is always the new branch connected to o"min going up to Amax- 

• As the temperature approaches its maximum value (T — )■ 1) the thin throat branch 
diminishes, and the part of the curve that remains coincides increasingly with that of 
the extremal one. 

• Because of the appearance of a maximum Amax and minimum Amin for temperatures 
in the range < T < Tf,, there exist values of the brane separation A for which there 
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are three possible phases. We will compare these in more detail in Ref. |20] . 



5.3 Analytical results 

We now consider what we can say analytically about the A versus ctq behavior either in the 



small temperature regime or the large ctq regime using the analysis of Section 4.3 for the 



branch connected to the extremal configuratior and the results listed in App. [b} 

We begin by considering the A versus ctq behavior for small temperatures T ^ 1. We 
first consider the minimum Amin of the A(cro) curve (for small temperatures this is a global 



minimum). The small temperature expansion for z'{a) in (4.34) can be used to compute the 



small temperature expansion of the brane separation (5.1) in the corresponding wormhole 



configuration. The result is given in App. [B| from which we quote 

A = Ao + T^Ai + C'(r^) 



(5.4) 



Here the leading term Ag is the extremal result (2.22) (with k = 1) and Ai is given in (B.3). 



Using now (5.4) we find that the minimum is found at 



(To = l + xf^ + 0{T^) 

1 , n-lfnir , V6r(-i)r(|)r(f) 



3\/6 1024 23/4 ^7r5/2 



-0.103 



(5.5) 



Note that the leading result ao = 1 reproduces the one found in Section 2.2 The corresponding 
value of the minimum separation distance is 



A, 



2V2^r(|) 
r(!) 



2vrr (i) 



(5.6) 



showing that for small temperatures the minimum is slightly lower than for T = 0. 

We turn now to the local maximum Amax for small temperatures. Using the expansion 
(B.2) of A in powers of up to and including the T^^ term we find that the ctq value at 



which the local maximum occurs is 



ao = air2/3 + a2Tio/3 + 0(ri8/3) 



(5.7) 



with ai and 02 given numerically by ai ~ 0.693 and 02 — — 0.00243j^ Note that the leading 
scaling (Tq oc T^/3 [yl (5.7) for the position of the maximum is consistent with the fact that 



it should be in between o"mm) which goes like T to leading order, and the ctq value in (5.5) 



corresponding to minimum Amin which goes like to leading order. The local maximum 

bi 



corresponding to (5.7) is 



A, 



J'2/3 ' 



(5.^ 



where 61 ~ 3.28 as computed using the expansion (B.2) to order T . Note furthermore that 



the contributions to each term in the expansion of fio in (5.7) come from each term of the 



^''Corresponding analytical results can be obtained for the branch connected to the neutral 3-brane, using 
App. [K\ but these will not be discussed in this paper. 



^Tb order in the expansion | |B.2[ ) we find ai ~ 0.707 and 02 — 0.0908. Going to order T one instead 
finds ai ~ 0.693 and 02 ~ -0.00243. 
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expansion of A in (B.2 ). This is the reason why we expand A up to T^^ instead of considering 
only the first two terms {i.e. up to the order T^) as we did in the determination of the 
minimum. 



We can also study what happens for ctq large. Here we can use (4.37) and in an expansion 



for small temperature we can compute the first correction, yielding 



A 



2ao + 



(5.9) 



We see that for large ctq the leading behavior of A is linear in ctq as in the extremal case (see 
Eq. 



6 Conclusions 

In this paper we have proposed a new method for D-brane probes in thermal backgrounds, 
studying in particular the thermal generalization of the Blon solution. To address this con- 
struction we have used the recently developed blackfold approach [TTl [IHl [H]. While this 
method was originally conceived and applied in connection with the approximate analytic 
construction of novel black hole solutions of Einstein gravity and supergravity in five and 
more dimensions, the results of this paper illustrate that it has a far broader range of appli- 
cability. 

In particular, the applications to new stationary blackfold solutions considered so far 
have mainly focused on black holes with compact horizons, and many new possible horizon 
topologies have been found. However, the approach is perfectly suited as well to describe the 
bending of black branes into other types of geometries. Consequently, it is the appropriate 
starting point to describe what happens to Blons when we switch on the temperature. Since 
the latter corresponds in the supergravity picture to an extremal D3-F1 probe brane system 
curved in a fiat space background [25] it is natural to base the thermal generalization on a 
non-extremal D3-F1 probe brane system curved in hot fiat space. The equilibrium conditions 
for such a configuration can then be computed from the blackfold equations. Moreover, 
as we have shown in this paper, the latter equations are in fact the natural non-extremal 
generalization of the DBI EOMs, providing an alternate (heuristic) derivation of the firmly 
established blackfold method. 

One may wonder how the method employed in this paper relates to previous works in 
which DBI configurations at finite temperature were considered. The commonly used method 
is to use the DBI action in a thermal background, e.g. hot fiat space or a Euclidean black hole 
background. However, this ignores the fact that, as soon as one switches on the temperature, 
the probe itself should be replaced by a thermal object. We believe the method developed in 
this paper is able to correctly describe this, as argued in more detail in the Sec. |3.3[ One can 
speculate that this new perspective on finite temperature D-brane probes might resolve the 
discrepancies between gauge theory and gravity results found in [U [15] . 

Another point worth emphasizing is that we have presented a thermodynamic action 
(Gibbs free energy) which gives the extrinsic blackfold equations describing the thermal gen- 
eralization of the Blon. This action reduces to the DBI Hamiltonian in the zero temperature 
limit, and may hence be regarded as the finite temperature/closed string analogue of the 
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corresponding DBI Hamiltonian. We also note that, as seen more generally in the blackfold 
approach, the thermodynamic origin of the action implies that the (mechanical) extrinsic 
blackfold EOMs are equivalent to requiring the first law of thermodynamics. The latter can 
be also seen as entropy maximization for given mass, so that this may be considered as a con- 
crete manifestation of an entropic principle governing equations that originate from gravity 
(see [26]). 



From the action we were able to obtain the explicit solution (4.26) for the slope of the 
embedding function describing the brane profile of a thermal D3-brane configuration with 
electric flux ending in a throat. We showed that there are two branches of solutions, one 
connected to the extremal configuration and the other connected to the neutral black 3-brane 
configuration. In most of our analysis we focused on the former branch. We have discussed 
the resulting finite temperature wormhole configuration of N D3-branes and parallel anti-D3- 
branes connected by a wormhole with F-string charge. We found that the finite temperature 
system behaves qualitatively different than its zero-temperature counterpart. In particular, 
for a given separation between the D-branes and anti-D-branes, while at zero temperature 
there are two phases, at finite temperature there are either one or three phases available. 

Moreover, from our results in Sec. [5] it seems that for small temperature and large enough 
(To the non-extremal Blon solution is well- approximated by the extremal Blon solution. We 
take this to mean that in this range using the (abelianized C/(l)^) DBI action as a probe of 
hot flat space could be a good approximation to our thermal D-brane probe. This would be 
interesting to examine further. It also illustrates that there are certain regimes where the two 
methods give different results which means that using our new method can change the results 
both quantitatively and qualitatively for certain regimes (e.g. the "thin throat" branch of the 
extremal solution becomes two branches with a maximum A). 

We left two important subjects for further considerations in the forthcoming paper |20j . 
The first subject is about the thermodynamics of the three branches of solutions that we found 
in Section [5} This is done by comparing the free energy for the branches in the canonical 
ensemble. The other subject is the apparent non-existence of a thermal generalization of 
the infinite spike solution as also seen in Section [5] For this we consider a different type of 
generalization by matching the supergravity solution of k non-extremal strings to our thermal 
D3-brane configuration with electric flux ending in a throat. 

In addition there are various open problems that would be interesting to pursue, on which 
we briefly comment. 

First of all, it would be interesting to generalize the solution of this paper to a thermal 
AdS background or AdS black hole backgrouncj^ This could have potential interesting appli- 
cations in the AdS / CFT correspondence and shed light on the above-mentioned discrepancies 
|14[ [T5] . More generally, since our method entails a new approach for D-brane probes in ther- 
mal backgrounds, it would be interesting to revisit other previously studied cases in which 
the classical DBI action is used for the D-brane probe in the thermal background. 

We have focused in this paper on the thermal generalization for the Blon in the case 
of D3-branes with electric flux. The construction is readily generalized for Dp-branes by 
starting with the EM tensor of the non-extremal Dp-Fl brane system. It would be useful to 



See Refs. [271 128] for applications of tlie blackfold approach in AdS backgrounds. 
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investigate this more general case, which could perhaps involve some qualitatively different 
features depending on p. 

The non-extremal D3-F1 brane system was treated in our construction at the probe level, 
ignoring backreaction effects. It would be worthwhile to try to go beyond the probe approx- 
imation and include such effects in a perturbative expansion. A scheme for this, based on 
matched asymptotic expansion, has been developed within the blackfold approach and al- 
ready been successfully applied to specific classes of black objects [171 IHl HI] • Furthermore, 
the blackfold approach has also given a powerful tool to study time evolution and stability 
\17\ [29] . Considering those for the thermal generalization of the Blon would be interesting as 
well. 

A more ambitious open direction would be to find a first principles derivation from string 
theory of the action describing thermal D-brane probes, generalizing the DBI action deriva- 
tions of [Ij. This would already be interesting to attempt to do for the case studied in this 
paper, for which we have found a thermodynamic action. We also point out that the analysis 
of this paper has been at tree-level in string theory. It might be possible to examine how the 
conclusions are affected by including one (or higher)-loop effects. 

Finally, we note that our method is not confined to D-branes but can be used more 
generally for all types of brane probes in thermal backgrounds, for example one could study 
M-brane probes in M-theory or NS5-brane probes in string theory. 
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A Analysis of branch connected to neutral configuration 



In this appendix we discuss some relevant limits of the solution branch (4.23) which is con- 
nected to the neutral 3-brane. 

We begin by considering the solution for small temperature. In this case, we find that in 
the limit of small temperatures we have 

A 1 

cosh^ a = 1 + — cos^ 6+— cos^ 6 + 0(cos^ 6) (A.l) 
27 243 

so that after some algebra 

--» = ^==(l + A4f«) (A.2) 

with the next correction being of order T^^{1 + k^/uq)^. The corresponding brane separation 
for the wormhole configuration then becomes 

2 V7rr(5/4) A ^ 8 ^2 
T(374r"°l'^27a-- 



^= ^Zo)r -o{l + ^-,f'] (A.3) 



32 



It is also interesting to consider here what happens when the D3-brane charge goes to 
zero, i.e. we have a neutral 3-brane with = 0. In some sense, this can be viewed as the 
opposite limit of the extremal limit connected to the first branch, which we examined above. 



Prom (4.4) we see that as — t- we have that k — )■ oo, so that from (4.20) we find 

ji4 



COS 5ia) = „ 



where we have defined 



9^ fcTpiT^ 
16^/3 T2 



(A.4) 



(A.5) 



D3 



Note that in the final expression the A^-dependence has canceled out, after substituting the 



definitions of k in (|4.4|) and T in (|4.17|). From (|A.4|) we immediately read off the lower bound 

(A.6) 



on a (and hence o"o) for a given temperature 



Substituting now (A.4) in the the relevant solution (4.23) for cosh a and using this in (4.25) 
we obtain 



-z'ia) 



1 + 



8 



(A.7) 



(A. 



Prom this we then find the corresponding brane separation as 

V^r(5/4) A ^ 8 f 8 
B A for small temperatures 



Here we show explicitly the small temperature expansion of A, i.e. the separation between 
the two systems of D3-branes. Using 



along with (4.11) and (4.22) we can expand A for small T up to the order T^^ 

A = Ao + T^Ai + f ^Aa + T^^Ag + 0{f^^) , 



where Aq is given in equation (2.22) and Ai, A2, A3 are computed to be 

Ai _2yfr(|) 



+ 1) 



-yVy^ + l2Fi 



1 5 1 
'2'4'"4' 



(B.l) 



(B.2) 



(B.3) 



A_2 



648yi3r(-i) 
(/ + 1) (4/ + 7) 2F1 



2 (2/ + 5/) 2F1 



1 5 
'2' 4' 



11 5 

'2'4'~4' 



31 (y^^ - 6/ + 11/ + 10) 



(B.4) 
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A3 _ x/fr(i) {y' + l) 

spk 22680yi9r(-|) 

> 3/2 



[7 (-42^^*^ - 2/^2 + 65/ + 145/ + 105) 



6/ (/ + 1)'^' (I71yi2 + 326/ + 294/ + 385) ai^i 



119 1 



+ 3/ (/ + 1)^/^ (342yi2 + 386/ + 504/ + 385) ^Fx 



2'4' 4' / 
15 9 1 



(B.5) 



2'4' 4' / 



where we defined y = o^j ^/k and 2-^1(0, c; z) is the hypergeometric function. 
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